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Abstract 

The density of states for the Schrodinger equation with a Gaussian random potential 
is calculated in a space of dimension d = 4 — e in the entire energy range including the 
vicinity of a mobility edge. Leading terms in 1/e are taken into account for N ~ 1 
(N is an order of perturbation theory) while all powers of 1/e are essential for N ^> 1 
with calculation of the expansion coefficients in the leading order in N. 

PACS numbers 03.65.-w, 11.10.HL 71.23.An 

It is a common belief!! that the average density of states has no singularity in the 
point of an Anderson transition in contrast to the conductivity and localization length. 
Nevertheless its calculation is of principal interest because all known methods fail in the 
vicinity of the transition. Another reason is that the density of states and the conductivity 
which are determined by the average Green function (G(x,x')) and the correlator (G R G A ) 
correspondingly are not completely independent quantities. In "parquet" approximation 
difficulties in both cases are of the same character and are connected with the problem of 
"ghost" pole!. On the other hand, to satisfy the Ward identity it is necessary to provide 
a strict correspondence of diagrams taken into account in the calculation of self energy 
and the irreducible vertex in the Bethe-Salpeter equation!. Thus no approximation for 
conductivity can be self-consistent before the corresponding approximation for density of 
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states is formulated. This problem was ignored in all existing theories! except the symmetry 
approach suggested in Ref. 5. 

For weak disorder the mobility edge is displaced near the bare band edge and the random 
potential can be regarded as Gaussian because the averaging is possible on scales smaller 
than wavelength and larger than the distance between scatterers (the so called Gaussian 
range of spectrum!) . The calculation of the average Green function for the Schrodinger 
equation with a Gaussian random potential is reduced to the problem of a second-order 
transition with n-component order parameter (p = (ifi,ip 2 , ■■■,fn) in the limit n — > 001. 
Then the coefficients in the Ginzburg- Landau Hamiltonian 

H{tp} = J d d x (ic|V^| 2 + \kM* + \uW\^ (1) 

are related with the parameters of the disordered system by the relations 

c = l/2m, Kq = -E, u = -a*W 2 /2, (2) 

where d is the dimension of space, m and E are the particle mass and energy, do is the 
lattice constant, and W is the amplitude of the random potential (we set c = 1, = 1 in 
what follows). The "wrong" sign of the coefficient of \(p\ A leads to unreacheability of Wilson 
fixed point in renormgroup equations! and to the problem of "ghost" pole in "parquet" 
approximation!. So the possibility of (4 — e)-expansion was in question for many years! and 
the number of suggestions on a value of upper critical dimension different from four was 
proposed0~lll. The progress was achieved in recent author's papers&ill, where the proper 
treatment of factorial divergency of perturbation series was shown to be necessary. Here we 
report the results for (4 — e)-dimensional case with the details of calculation to be published 
elsewhere^. 

In a four-dimensional space the structure of the perturbation series for the self-energy 
E(p, k) at p = has a form0 

00 N / A\ K 
£(0,k)-£(0,0) = K 2 5>^£^ In- , (3) 

JV=1 K=0 ^ K ' 



where k is the renormalized value of kq and A is the large-momentum cutoff parameter. The 
analogous expansion for d = 4 — e has the form 

k 2 + E(0, k) - S(0, 0) = k 2 F («) = f; (uA-y J2 A %(e) [ (A/k)£ ~ 1 

where A^(e) are the regular functions of e, 



(4) 



^(c) = E ^' L e L (5) 

L=0 

and A(j(e) = 1. The expansion (4) takes account of the fact that Y is a homogenious poly- 
nomial of degree N in A _e and K~ e , as follows from the dimensional analysis, and that the 
expression (4) should be reduced to (3) in the e — > limit. 

The quantity Y satisfies the Callan - Symanzik equationEl, which follows from its relation 
with the vertex r^'^El: 

+ W(g ,e)— + V(g ,e))Y = 0, (6) 



\<91nA dg 

where go = uA~ e , V(go, e) = ^(go,e), and the functions T^' 2 \ W(go,e), and i] 2 (go,e) are 
defined in Ref. 18. Introducing the expansions 
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W(g ,e)= E W M (e)g M = E E Wmm'9o^\ 

M=l M=l M'=0 (n\ 

oo ,,0000 \ / 

^(</o,<0 = E V^(e)(^= E E VW,M^o M e M 

M=l M=l M'=0 

and substituting expression (4) into Eq. (6), we obtain a system of equations for the functions 

Am ■■ 

N-K 

(K + l)A«+\e) — (N — K)eA K N {e) - £ [{N - M)W M+1 (e) + V M (e)]A N _ M (e) (8) 

M=l 

and for the coefficients A^ ,L : 

N-K L 

(K+l)A* +1 > L = (N-K)A«' L -\l-5 L>0 )-Yl E [(N-M)W m+ i,m'+V m , m/ }A^ i m ' . (9) 

M=l M'=0 

In the standard procedure of e-expansionl a few first terms in series (5) are retained. 
In M-th order in e the coefficients A^~ K,L with K + L < M — 1 are necessary for which 



the closed system of the difference equations is followed from Eq. 9. Initial conditions to 
this system and coefficients W 2t o, Vi,o, ■ ■ ■ can be derived if a few first orders of perturbation 
theory are calculated. By separating out the leading asymptotic term in N, it is easily 
proved by induction that 

A N N - K ' L = C« +L A»- K - L , (10) 

lJV-iC _ I tj/ \N 



h TjN-Po) (-W 3 , ) K (N\nN) K 

v 1 2 ' 0J r(iv + i)r(-/3 )(-^ 2 ,o) 2X k\ : 



where (3q = — V^o/W^o, and the value of the first few coefficients in the expansion (7) equil0 

W 1 (e) = -e, W 2j0 = K A (n + 8), W 3}0 = -3K 2 A {3n + 14), V 1>0 = -K A (n + 2). (11) 

with Ki defined in Eq. 13. For the parquet coefficients the result (10) is exact. The 
coefficients (10) do not possess the factorial growth which for u < is responsible for 
nonperturbative contribution resulting in existence of the fluctuational tail of the density of 
stateSH. This is a reason why the usual Wilson method does not work for u < 0. 

The approximation giving asymptotically exact results for small e is as follows: In the 
lowest orders of the perturbation theory it is sufficient to retain in expansion (4) only leading 
order in 1/e; for large iV the lowest powers of 1/e should be taken into account, since the 
corresponding terms grow rapidly as N — > oo, but the leading order in N is sufficient for the 
expansion coefficients in Eq. 4. Information about the coefficients Ajy(e) for large N can be 
obtained by the Lipatov method00. The iVth order contibution to S(p, k) is calculated in 
a close analogy with the d = 4 caseEl and has a form 

[S(p, k)] n = c 2 u N T (N + b)a N J d\nR 2 R-\4,l) Rp {4,l)-Rp- 

o 

•exp (-Nf(KR) + NelnR + 2K d h{KR) 1 - ) , (12) 



where 



-3K A , b=^, C2 = c(3K A ) 7/2 
1 . x 



f(x) = ~\{C + 2 + lnvr) - 3x 2 (c + - + In |) , (<j> c ) 3 p = 8 ■ 2 1 /\ 2 pK 1 (p), (13) 
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h(x) = 7 4 exp(/(x)), h = yS 4 , S d = 2n d / 2 /T(d/2), K d = S d (2ir)- d , 

C is Euler's constant, Ki(x) is a modified Bessel function, and c is a constant of order unity 
defined in Eq. 114 of Ref. 16. Representing the result (12) in the form of expansion (4), we 
have 

A«(e) = 5 2 r (N + b)a N C§Jd\nR 2 R- 2 (e + ^± e m)-^R 



K 



■exp \-Nf{R) + Ne\nR + 2K d I 4 (R) 1 ^ j , (14) 

where c 2 = c 2 ((f) 3 c )l « 3.44 • 10~ 2 for n = 0. 

As discussed in Ref. 16, for d — 4 the Lipatov method reproduces the coefficients 
in (3) well only for K <C A, since they decrease rapidly with increasing K and accuracy 
~ 1/A of the leading asyptotic expression is limited. Analogous phenomenon takes place in 
(4 — e) dimensions: Formula (14) is valid for all K if Ae 3> 1 and only for K <C A if Ae < 1; 
under these conditions the coefficients (14) satisfy Eq. 8, where only the term with M = 1 
is retained in the sum, which is possible for large values of A in view of factorial growth of 
A§(e). Since the system of equations (8) determines A§(e) for K > if A° N (e) are given, 
and the latter are well reproduced by the Lipatov method, one can determine A^(e) in the 
region 1 <C A < 1/e for all K. 

Retaining in Eq. 8 only the terms with M = 1 and M = 2 and introducing the quantity 
Xn,m by definition 

4?M = W^ fjFTW - 2- (15) 

where /3(e) = — Vi(e)/W2(e), we have an equation 

Xn,m = (Jm + <5m) Ajv,m+i , (16) 



where 



l M = h M + e-* , 5 M =^e~ 2i P (17) 
hM --W 2 (e) A° M (e) M-fl(e) ' ^ ~ W 2 (e) A%(e) P( 6 )) (18) 



and is the displacement operator on the distance —1 affecting the both arguments. 
The equation (16) can be formally solved by straigtforward iterations taking into account 
the boundary condition X^n = 1. The result can be expanded in powers of 5m, terms 
containing a few number of operators 5m being essential. These terms can be explicitly 
calculated by induction. For M ~ 1 we have a result 

A N N ~ M {e) = -L e N - M c 2 T{N-{3{e))a N (t/2n) 1/2 e f~(™mN-i) + i/t. 

■ I™ dxe-^-^-^x^+P-^Jix ) , (19) 
Jo 



where 



J(N) = J d In R 2 R~ 2 exp f —Nf(R) + NelnR + 2K d h{R) 1 ^ j , (20) 



and 



t _ ea 3e W 3 (e) e -+o 3n + 14 

W 2 (e) + 8 ' /o °"alf 2 (e) * n + 8 1 J 

which is valid in the region Nt > 1 or 1 - Nt < e 1/2 . For < 1 and iV - K < IniV we 
have a result of type (5) with A^ L given by (10). Other regions of parameters can be also 
investigated but they do not make essential contributions to the sum of the perturbation 
series. 

Two contributions are important in the sum of Eq. 4: (a) nonperturbative contribution 

P(0, n)} nonpert = ir (K 2 ) = m~c 2 K 2 {^/au) b e - K '/ au F (n'/au) (22) 

(the limit A — > oo is taken) arising from the region of large A" and obtained from Eq. 4 by 
summation over A" from an arbitrary finite Nq to infinity, if the coefficients A^(e) given by 
(19) are written in the form 

A#(e) = 5 2 r(iV + b)e N a N F(N) (23) 

(the formula (46) of Ref. 16 is used); (b) the quasiparquet contribution arising from the 
terms with coefficients A^~ K ' L with K ~ L ~ 1 given by (10) 



quasiparq 



A -| —uK~ e In A 

W 2 , 



A) 



K~ e - A" 



c 



A ee 1 + W 2 , u . (24) 



To logarithmic accuracy, the quantity A in the logarithm can be replaced by its minimum 
value A ~ elne (determined by Eqs. (26) -(30) presented below), and in the limit A — > oo 
the result (24) can be written in the form 



[Y(K) guasiparg = [1 + W 2 , uk 7e] A 

which differs from the parquet result^ only in that u is replaced by u. It can be proved^ 
that such replacement occurs in all parquet formulas employed in calculating the density of 
states^. 

The rest of the calculations are similar to those described in Ref. 15. The damping r, 
the renormalized energy E, and the density of states v are determined in parametric form 
as functions of the bare energy E B by the equations 



u = u 



W 3 ,o 

W? 



eln A 



2.0 



(25) 



r = r c [i + —) sin^, E = -r c [i + —j cos^, 



(26) 



v = — — 1 + — 

4vr u V 2 



2\~ 1/4 / V 
1 + — sm [<p + — 
ex J V Ax 



1 + — smU+f 
ex / V 4x 



w(l + 2ex) . / if 



Rq 



2(1 + ex/2) 



r, 



'8iG|«| 



2/e 



d d k 1 



1/2 



E c ~2u, 

' (27i) d k 2 



Ro 



(27) 



(28) 



(29) 



e ) J ['Att)"' K" \3\n(l/e) J 

where x(ip) is a single- valued function in the interval < if < tt, similar to the function 
shown in Fig. 2 of Ref. 15 and determined by equation 

- 4x '/ 3 <p(l + 2ex) 



x 



1/4 



-/(x) cos ■ 



(30) 



where 
I 



/5\ 1/4 /7r y- x 1/2 rv 
(#) = C 2 (-) ( — ) e -/oo+/oo(l+^/2)ln[A(l+ex/2)/t] y § (ex/2t-zf ^+/3-/oc (l+ex/2) 

(31) 
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The expressions (26-30) are simplified drastically in two overlapping regions. For large \E\, 
when x ^> In (1/e), the right-hand side of (29) is small and the quantity (p is close to or ir, 
so for r(E) and v(E) we have asymptotic expressions 



r(E) 



-eE 



(E/r c y/* - i 



-i 



r (E) i-(\E\/r c ) 



-6/2 



-1/4 



£ > r 



(32) 



u(E) 



2 



r (E) 



i 



Rk l\E\ 



-6/2 



1 



\E\ 



-6/2' 



4vr|w| \ 2 \ r c J 

(r (E) = Jo(|/c| 2 )) which give an illusion of ghost pol 
Eq. 29 takes the form 



1/2N 



(33) 



For small \E\, when x 5; e -1 / 2 



sin (93 + ip/Ax) = I(0)cos(<^/3) — j- 



-4x/3 



x x/4 



J(0) 



/ 1x17/12 

7 / 12 (\n-\ 



e) 



(34) 



and the ghost pole is shifted from the real axis into the complex plane on the distance 
~ eln(l/e). 

For large positive E Eq. 32 gives the inverse relaxation time appearing in the kinetic 
equation while for large negative E the damping r becomes purely nonperturbative. The 
function u[E) for large positive E goes over to the density of states of an ideal system, and 
for large negative E the following result is obtained for the fluctuational tail 

1 \ 1/2 



u(E) = ^r (E)\Er/Hn-L = ~c 2 K 4 f^m-M n~ : ''\E\ { ''~ 2 "" 2 
71 Rq V 3 Rq. 



■ exp 



h 


E\ 


6/2 


4 


u 





(d+l)/2 



h(Ro) 


\E 


6/2 \ 


4 


u 


Rq J 



(35) 



the energy dependence of which is identical to that obtained in Ref. 22 and 23 and corre- 
sponds to Lifshitz's lawi; the divergence in the limit e — ► is removed for a finite cutoff 
parameter A. 

It is interesting, that for ex <C 1 formulas (26-30) have the same functional form as in 
the d = 4 case@, i. e. the behavior of all physical quantities in the vicinity of a mobility edge 
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is effectively four-dimensional. As in Refs. 15 and 16, the phase transition point shifts into 
the complex plain, and the density of states has no singularities for real E in accordance 
with widely accepted but not proved ideas. 
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